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Abstract 



We analyse a new class of statistical systems, which simulate different sys- 
tems of random surfaces on a lattice. Geometrical hierarchy of the energy 
functionals on which these theories are based produces corresponding hierar- 
chy of the surface dynamics and of the phase transitions. We specially consider 
3D gonihedric system and have found that it is equivalent to the propagation 
of almost free 2D Ising fermions. We construct dual statistical system with 
new matchbox spin variable G^, high temperature expansion of which equally 
well describe these surfaces. 



1 . The correspondence between spin configurations and the surfaces of interface 
allows to define different theories of random surfaces on a lattice. Ising ferromagnet 
is equivalent to a system of random walks with the energy which is proportional to 
the length of the path on the two-dimensional lattice and to the area of the surface 
in three-dimensions § % §, |, §, § % 



In the recent articles [|K| [L3| the authors formulated a new class of statistical 
systems, whose interface energy is associated with the edges of the interface. Spin 
systems introduced in |T0| , [I3|| represent a system of random walks, but with the 
energy which is proportional to a total curvature of the path on the two-dimensional 
lattice and to the size of the surface in three dimensions fT3| , [T5| , [UJ . One motivation 
for the study of those statistical systems is to have well defined physical system of 
random surfaces which is relevant for the description of the gauge field dynamics 
and of the QCD 0, |15|, |16| . 

In order to study the dynamics of the surfaces with linear-gonihedric action 
[O , |15| , |16| we analyse statistical systems with topological and area action and 



compare statistical properties of these systems from the viewpoint of surface dy- 
namics. The energy functionals on which these theories are based represent a real- 
ization of the hierarchy of geometrical functionals found by Steiner and Minkowski 

The imposed structure has an important influence 



[0, [Z§ fSJ [27||. 

on the physical properties of the system and produces corresponding hierarchy of 
the surface dynamics and of the phase transitions. We have found that 3D goni- 
hedric system is equivalent to the propagation of string with almost free 2D Ising 
fermions. We construct dual statistical system with new matchbox spin variable G^, 
high temperature expansion of which equally well describe these surfaces. 

2 . The origin of the interaction hierarchy lies in Steiner's geometrical idea about 
parallel manifold M p [19, [14], 21 , ^] . His formula represents the area of the parallel 



surface M p , situated on the distance p from M, in terms of the area S, mean integral 
width A and the topology T of the initial convex surface M 



S p = I RxRzdu + p ( (Ri + R 2 )duj + p 2 [ dco, 
Jm Jm Jm 



(1) 



where -Ri^ are the principal curvatures of the surface M and du is the area element 
of the spherical map. All this functionals, S(M), A(M) and T(M) properly defined 

mi 



for nonconvex case 



mm 



produce the systems of interacting surfaces with 
a very deep hierarchical structure. 

At the first level one can define topological-gonimetric systems on the lattice, 
these are ID Ising ferromagnet, 2D gonimetric walks and 3D gonimetric surfaces 
PL 1137 



H lFing = ~ a(T '' H gonimetric = ~ H gonimetric = ~ Yl OOOOOOO^ 

links plaquettes boxes 

(2) 

where the spin variables are attached to verteces of the lattice. Using high temper- 
ature expansion one can see that all topological systems (2) have the same partition 
function and are in a completely disordered phase |TB| 



~Pf gonimetric ~ lfl(u + UJ 1 ), (3) 



where uj = e' 13 . 

Increasing the dimension of the lattices by one and leaving the Hamiltonians 
(2) without changes, one can see that the system (2) "moves" to next level of the 
hierarchy and describes now the system with gonihedric action A(M) = J2<ij> \Xi — 
Xj\ ■ 1 7T — ctij\ which is proportional to the linear size of the system |T4|, |T5], [TH ]. 



There is mutual connection between gonimetric and gonihedric systems. Indeed 
in the recent article ]28| the authors have found that the action A(M) has an equiv- 
alent representation in terms of total curvature k(E) of the polygons which appear 
in the intersection of the two-dimensional plane E with the given surface M . This 
curvature k(E) should then be integrated over all planes E intersecting the surface 
M. This result directly connects topological-gonimetric system in two dimensions 
with gonihedric system in three dimensions. This connection is universal for any 
dimensions. 

Indeed the intersection of the two-dimensional plane E with the polyhedral sur- 
face M is the union of the polygons 

P 1 (E),...,P k (E) (4) 



and as it was shown in |[28|| , the absolute total curvature k(E) of the polygons (4) 
in the intersection is 



fc(JE0 = £ k-agl, (5) 

i=l <i,j> 

where afj are the angles of the polygons (4) and they are defined as the angles in 
the intersection of the two-dimensional plane E with the edge < Xi,Xj >, that is 
with the dihedral angle . Particularly for the plane E which is perpendicular 

to the edge < Xi,Xj > we will have that afj = azij. The meaning of (5) is that it 
measures the total revolution of the tangent vectors to polygons (4) Pi(E), Pk{E) 



2S, 30, 31, 32]. The summation in (5) can be extended to all edges of the surface 
M if we take afj = ix for the edges which are not intersected by the given plane 
E. Finally to get an action A(M) one should integrate the total curvature k(E) (5) 
over all intersecting planes E 



MM) = TJ {E) k(E)dE - (6) 

where the integral is extended over all planes{i?} intersecting the surface M. The 
measure dE in (6) is equal to 

dE = dpdVt (7) 

where we define the plane E by its normal vector n = (sin 9 cos 0, sin 9 sin 0, cos <$) 
and its distance to fixed origin. dQ denotes the area element of the unit sphere 
corresponding to the end point of the normal vector ft. In the next section we will 
apply (6) to lattice surfaces. 

3. Singular surfaces {M} of interface Q on a cubic lattice Z 3 can be consid- 
ered as a collection of plaquettes with the restriction that only an even number of 
plaquettes can intersect at a given edge ( 2r = 0,2,4.) and that only one plaquette 
is on a given place [S [TU], [LlJ . 



As a first step we should find the same representation (6) for the action A(M) 
on a cubic lattice Z 3 . This can be easily done if we introduce a set of planes 
{E x }, {E y }, {E z } on the dual lattice Z\ which are perpendicular to x,y and z axis 
correspondingly. 

These planes intersect a given surface M and on each of these planes we will 
have an image of the surface M. Every such image is represented as a collection of 
closed polygons {P(E)} appearing in the intersection of the plane with surface. 

The energy of the surface M is equal now to a sum of the total curvatures k(E) 
of all these polygons on different planes 



A(M)= £ k(E). (8) 

over all dual planes{E x ,E y ,E z } 

On the lattice the total curvature k(E) (5) is simply the total number of polygons 
right angles. In general, the surface on the lattice may have self-intersections and one 
should associate the energy with self-intersection edges of the surface and therefore 
with the self-intersections of polygons. Depending on how we will count the right 



angles in the self-intersection vertices we will get different theories Rig , |10[ y|. 

In the present article we will consider the case when the self-intersection coupling 
constant k is zero [[DJ. This means that when we compute the total curvature k(E) 
of the polygons in (8) we should ignore the right angles at the self-intersection points. 

In terms of Ising spin variables ay the Hamiltonian of this system has the form 
il 



Hgonihedric ~ (Tp<7f + ^0 '^ +s+ pO (2a) 

r,a,/3 

—* 

where r is a three-dimensional vector whose components are integer and a ,(3 are 
unit vectors parallel to axis. The gonihedric system (2a) has an extra symmetry: 
one can independently flip spins on any combination of planes (spin layers). 
4. With (8) the partition function of the system can be written in the form 

Z((3)= £ exp{-(3j2k(E)}. (9) 

over all {M} {E} 

Let us represent the sum in the exponent as a product 



expi-^KE)} = J] = J] e _/JW II ^ KEv) II ( 10 ) 

{E} {E} {E z } {E y } {E x } 

It is true that these products are not independent for a given surface M and our 
goal is to express the initial energy (8) and the product (10) in terms of independent 
quantities, that is only through the product over all two-dimensional planes in one 
fixed direction, let's say through {E z }. 

The question is: what kind of information do we need to know on these planes 
{E z } to recover the values of the total curvature k(E x ) and k(E y ) on the planes 
{E y } and {E x } for the given surface M ? 



With this aim let us consider the sequence of two planes E\ and E z +1 and denote 
by Pi the polygon- image of the surface M on the plane E\ and by P i+ i the polygon- 
image of M on the plane E l z +1 . 

With this one can compute the contribution to the total curvature k(E x ) which 
comes from the pieces of the polygon- images on the planes {E x } which lie between 
two planes E\ and E z +1 . This contribution is equal to the number of bonds of the 
polygons Pi and P i+ i which are parallel to the x axis without common bonds. 

In the same way the contribution to the k(E y ) from the polygons on the planes 
{E y } which lie between E l z and E z +1 is equal to the length of the non-common bonds 
of the Pi and Pi+i which are parallel the y axis. 

Therefore the total contribution to the curvature k(E x ) + k(E y ) of the polygons 
which are on the perpendicular planes between E\ and E l z +1 is equal to the length 
of the polygons Pi and P i+ i without length of the common bonds 

i(p t ) + i(p +1 )-2-i(pnp +1 ). (li) 

This formula represents an important fact that the curvature of the polygons which 
lie on the perpendicular planes {E y } and {E x } is equal to the length of the polygons 
on the {E z } planes 

Contribution to k{E x ) + k(E y ) \ 

from planes {E x }, {E y } = l(P t ) + l(P i+1 ) - 2 • l(P n P i+1 ). (12) 
between E\ and E l z +1 J 

Where 1(1) + 1(2) - 2 • 1(1 n 2) = 1(1 U 2 \ 1 n 2). Now the energy functional (8) is 
reduced to an independent sum over the polygon loops only on the E z planes 

A(M) = J2 HE) = J2 k(Pi) + KPi) + KPi+i) - 2 • l(Pi n P i+ i) (8a) 

m {e z } 

which we will consider as a sum of the free action proportional to the total curvature 
plus the lenght 

A (P) = k(P)/2 + l(P) (8b) 

and of the interaction term which is proportional to the length of overlapping of 
the polygon loops Pi and P i+1 

A int (P, P i+l ) = -2 ■ l(P t n P i+1 ). (8c) 
With this formula one can represent the product (10) in the form 

Yl e -^) = Yl exp{-(3(A (P l ) + A (P l+1 ) + A int (P h P t+1 )} (13) 

m {e z } 

and the partition function as 



Z (P)= E Y[ex P {-P(A Q (Pi)+A Q (P i+1 )+Ai nt (Pi,P i+1 )}, (14) 

{Pi} on {El} i 



where an independent summation is extended over all closed polygons on the differ- 
ent planes E z . 

If we define transition amplitude from configuration Pj on the plane E\ to the 
configuration P i+1 on the plane E l z +1 as 

K{Pi, P i+1 ) = exp{-P(A Q (P t ) + A {P i+1 ) + A int {Pi, P i+1 )}, (15) 

then 

Z(P) = Y, II K(P h P i+1 ) (16) 

{Pi} i 

which we can interpret as the propagation of the polygon-loop in the z-direction. 
To compare with 3D Ising system |3B], [55], |3(| [J7|, M we will represent the 



energy of the interface in the form 

S(M) = W l (E) = E KPi) + <Pi) + <Pi+x) - 2 ■ s(P t n P i+1 ) (17) 

m {e z } 

with the free action A which is proportional to the total lenght plus oriented area 

A = l(P)/2 + s(P) (18) 
and interaction which is now proportional to overlapping area of the polygon loops 

A int = -2 ■ siPi n P i+l ) (19) 

where s denotes the area. We see, that compared with (8), the interaction is much 
more stronger. 

Description of the spin systems in terms of interface factors part of the symmetry: 
two spin configurations connected by the global Z 2 transformation correspond to the 
same surface. Description in terms of polygon loops (8a,b,c) factors the number of 
surface configurations by the factor 2^, when n = 0. This fact is connected with 
the layer symmetry of the system when k = (see section 3). 

The 3D Ising ferromagnet does not have this symmetry, therefore in (17), (18) 
and (19) we have to ascribe consistent orientation to every polygon loop on the E l z 
planes and then to sum over 2 N different orientations. 

5. Using (8a,b,c) let us define 

A{1, 2) = A (l) + A (2) + A mt (l, 2) (20) 
where Pi = % so that (15) becomes equal to 

if (1,2) = e~ p A{1 ' 2) . (21) 

Let us consider intermediate summation over all polygon configurations between El 
and El 

£ K(l, 2)K(2, 3) = £ e- p A(2 ' 3) (22) 

{2} {2} 



and represent the result in the form 



e -0A(i,3) . F (i )3 ) ) (22a) 

where 

F(l,3) = E e_2 ' /3 A ° (2) ^(1,2,3), (23) 

{2} 

and 

V(l 2 3) = e' 13 (~ A ^t(i,s)+A int (i,2)+A int (2,3)) ^4) 

The last term in (24) represents the interaction which is proportional to the over- 
lapping lenght. 

6. To simplify the system and to have crude approximation we will ignore the 
interaction term (24). In that case F(l,3) does not depend on the loops 1 and 3. 
Therefore denoting this amplitude simply by F we will have 

F = J2 e' w ' Ao{p) = Y, e~ mp)+2l{p)) (25) 
{P} {P} 

and that the partition function reduces to two-dimensional partition function 

Z (P) = F N . (26) 

So in this approximation the partition function coincides with the two-dimensional 
model (25) with the energy which is proportional to the sum of the lengths of 
the loops plus the total curvature. This system coincides with the sum of the 
2D Ising model Z{(3) = Y,{p} e_/3 ' an d of the 2D gonimetric walks [II], |T^] 
Z{(3) = D{p} e_/3 ' k<yP ^ ■ In short one can say that in this approximation the sys- 
tem is equivalent to 2D Ising model in which the paths are weighted by the total 
curvature. Our aim is to evaluate (14) in this approximation. 

T. The two-dimensional topological model (2), 2D gonimetric walk, has the 



partition function (3) and is in the disordered regime [13 |. 

For the three-dimensional case, we obtain an additional perimeter term to goni- 
metric walks (8b), (25) and this results to the change of the phase structure of the 
system in three dimensions. We can expect that because of the perimeter term, the 
linear system (6), (8) will show the phase transition in 3D which should be of the 
same nature as it is in the 2D Ising ferromagnet. 

To find partition function one has to represent the corresponding weights in 



terms of eight- vertex model || Bq, KJJ, 11]. The 2D Ising system has the weights 



UJ\ = 1, U 2 = W 4 , U 3 = kJ 4 = Co> 5 = U)q = OJj = UJ 8 = U> 2 , (27) 

where w = exp(—j3J) and = exp{—[3e^) and eg is the energy assigned to the £th 
type of vertex configuration (£j = 1, .., 8) e\ = 0, e 2 = 4J, e 3 = = e§ = e§ = ej = 



es = 2 J. The weights of the gonimetric walks with k = are [110, [13 



LO\ = L0 2 = ^3 = ^4 = 1, CU5 = Co>g = UJ7 = Co>g = UJ, (28) 



where uj = exp(—(3Q{ji/2)) and e\ = t<i = £3 = £4 = 0, e$ = e§ = e-r = eg = Q(tt/2). 
With this notations the system (25), (26) has the weights 



fix = 1, tt 2 



w 



w 



1^5 — fig — — fig 



• w 



(29) 



Following [J13, |6|, |(| [38|, 11| one can rewrite the partition function (25) as a product 
of the fermion operators and obtain 



-PfHonihedric = ^ M(l + W A ) 3 ~ 4«;V(1 - U 2 ) + Aw\l - U 2 ) COS 6 COS 

+2(w 6 + - 2wV)(cos0 + cos (j))]d9d(f>, (30) 

which confirms our expectation. One can also compute the partition function with 
the Kac and Ward || combinatorial approach. To each oriented closed polygon 
one can correspond a term of the matrix determinant and vice versa. For that we 
should define the matrix A((i,j)(X,Y)), where - indicates which points of 

the lattice are connected by the polygon loop and X = (Right, Left, Down, Up) = 
(R, L, D, U) - indicate the direction in which the bond joining the vertices of the 
loop is traversed and Y = (R, L, D,U) - subsequent direction of "motion". In our 
case RR = LL = UU = DD = w, RU = LD = UL = DR == wua, RD = 
LU = DL = UR = wua~ , where a = exp(— in/A) ||. The computation of the 
determinant can be easily done and coinside with the expression (30). We conclude 
from this that only boundary terms can change this result and that the system 
describes the propagation of almost free fermionic string. 

8. We have found also the dual system to (2a), high temperature expansion of 
which equally well describe those surfaces. The high temperature expansion of (2a) 
is 

Z((3) = J2 II ch(3 ■ {l + th/3 ■ (aaaa)}. (31) 

{<r} plaquetts 

Opening the brackets and summing over a one can see that only such terms produce 
nonzero contribution which contain an even number of plaquettes on every given 
vertex, therefore 

Z(0) = (2chi3) 3NS EW) S(E) , (32) 
{£} 

where the summation is extended over all surfaces {£} with an even number of 
plaquettes at any given vertex. The s(E) is the number of plaquettes of E, e.g. the 
area of the surface. 

Let us attach plaquette variables Up to each plaquette P of Z 3 

Up = -1 if P e M and U P = 1 if & M (33) 
The constraint on the plaquette variables Up in every vertex 



n u p =i > 

12 plaquettes incident to vertex 



(34) 



uniquely characterizes our set of surfaces {£}• Now one can introduce the group 
structure on this set of surfaces {£}. Let us consider two surfaces S 1 and S 2 and 
denote their plaquette variables as Up and Up respectively. Let us define group 
product of these two surfaces as 

U P = Up- Up. (35) 

According to this definition the set of surfaces {£} (34) forms an Abelian group 
G. The whole group G is a direct product of the local groups G^. This group 
G^ has four elements-elementary surfaces, G^ = {e(£), g x (i~), g v (£), w hh the 
multiplication table e • g X:VS = g x>vs ; g x -g x = g v • g v = g q ■ g q = 1; g x -g n = which 
follows from the fact that elementary surfaces are match box surfaces with different 
orientations and from the multiplication law (35). 

Any set of elementary surfaces e,g x ,g v ,g q distributed independently over the 
lattice Z 3 describes some allowed surface E and any given surface from {£} (34) 
can be uniquely decomposed into the product of G% 

£ = n-GV (36) 

This approach allows to describe the original surface £ in terms of a new independent 
variable G^ = {e(£), <7 X (6)> fl^CO) 9s{£)} which should be attached to center of the 
cube £ of the original lattice Z 3 or which is the same to the vertices £ of the dual 
lattice Z* 3 . 

The group G^ is an Abelian group of the fourth order and therefore has four one- 
dimensional irreducible representations E — {1, 1,1,1}, R x = {1,1, — 1, — 1}, R v = 
{1,-1,1,-1}, R' = {1, — 1, — 1,1} which express algebraically the "matchbox spin" 
variable G^. 

The dual Hamiltonian is nonhomogeneous in the directions Xi V an d £ 

where Xi V an d q are unit vectors in the corresponding directions of the dual lattice 
and 

H u+X = H{G f , G (+x ) = -J2*(fl • i2*(£ + X ), 

H u+n = H(Gs, G e+J? ) = -#>(£) • + rj), 

H^ + , = H(G i7 G c+? ) = -i? ? (0 • i? ? (£ + q). (38) 

The partition function of the dual system (37), (38) can be written in the form 
Z(f3*) = X){gj exp{— j3*Hduai}- One can check now, that the high temperature 
expansion of the dual system (37), (38) indeed coinsides with the low temperature 
expansion of the original system (2a) and provides a new realization of gonihedric 
system on the lattice. This dual representation can be compared with the dual 
representation of the 3D Ising ferromagnet in terms of 3D Gauge Wegner system 

!• 
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